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' The time evolution of a bounded quantum system is considered in the framework of the orthog- 

, onal, unitary and symplectic circular ensembles of random matrix theory. For an N dimensional 

^\ • Hilbert space we prove that in the large N limit the return amplitude to the initial state and 

the transition amplitude to any other state of Hilbert space are Gaussian distributed. We further 
compute the exact first and second moments of the distributions. The return and transition prob- 
• abilities turn out to be non self-averaging quantities with a Poisson distribution. Departures from 

| this universal behaviour are also discussed. 
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I. INTRODUCTION 



During the last decades much theoretical progress has been made concerning the physical properties of complex 
quantum systems. In the 1950s and 1960s, following the pioneering work of Wigner, a statistical theory was developedtl 
to interpret the basic features of highly excited nuclear states. Known today as the random matrix theory (RMT), it 
has had an impressive success in the description of a large class of systems, ranging from nuclei and atoms down to 
microwave chaotic cavities or transport properties of electrons in metallic or ballistic disordered quantum dots (for a 
review of the developments see for instance Refs. ||^,^). 

Most past studies focus on the spectral statistics of complex systems based on linear wave dynamics. Our purpose 
here is to concentrate on the time evolution of such systems. We are particularly interested in specific quantum (or 
wave) manifestations in the dynamics. If at t = the system is prepared in a given state |"0(O)), a quantity which 
has been extensively studied is the average of the return probability |C(t)| 2 = |(^i(0)|^(i))| 2 (sometimes called the 
^\ , survival probability). This function, which is accessible from experimental measurements as the modulus squared of 
the Fourier transform of the spectrum, was studied by Leviandier et alEI in their analysis of spectral correlations in 
molecular physics. Their basic finding is that the spectral rigidity is manifested as a "hole" in the average of |C(i)| 2 
for times shorter than the Heisenberg time (proportional to the inverse of the mean level spacing), this test being 
more robust than other ways of testing the spectral correlations. The technique was also applied to the nuclear data 
ensembled and the spectra of superconducting microwave cavitieso. From the theoretical point of view, a rigorous 
interpretation of this phenomenon in the context of RMT was given in Refs. SJlQ (see also DTI). The correlation 



. hole was also obtained using supersymmetry techniques and interpreted as a quantum dynamical echo in the time 
evolution of a wave packet within a disordered mesoscopic sampleEJ, and in acousticsEil In Ref. |l2| a direct dynamical 
experimental test of the correlation hole based on the transient currents through a disordered dot was also proposed, 
O ■ but such a direct test is still lacking. 

Our purpose here is to describe the statistical properties of the return probability as well as those of more fundamen- 
tal quantities such as the return amplitude C(t) and the transition amplitude to a different state \x), T(t) = (x\ip{t)}. 
Our basic assumption is that the system can be described, from a statistical point of view, by a random matrix 
belonging to one of the three well known symmetry classes of RMT. In this respect, and due to the Hilbert space 
rotational invariance of these ensembles, the theory applies to the relaxation process of a closed system for which the 
rate of exploration of phase space is sufficiently high in order to avoid localization effects. We show in the following 
that in the large- N limit the real and imaginary part of C(t) and T(t) are Gaussian distributed for any time. It then 
follows that their modulus squared (the return and transition probabilities) are Poisson distributed and that moreover 
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these are non self-averaging quantities (in the sense that their fluctuations are of the same order of magnitude as 
their mean value in the large- iV limit). This provides incidentally a proof of the results of Ref. [If] where the Poisson 
distribution was conjectured. Finally we obtain an exact expression, valid for any N, for the form factor and thus 
for the average return probability. This generalizes the asymptotic results obtained in Refs. p|, ^o|Jl^ . An intuitive 
interpretation of the results obtained in terms of a correlated random walk process in the complex plane is also given. 

Deviations from the universal behaviour associated to a RMT description are expected to occur in particular for 
short times. These effects, which are not initially included in our description, are briefly discussed in the last part of 
the paper. Our results are moreover obtained for the circular ensembles of RMT. In the-context of dynamical systems 
these ensembles are closely related to the quantized version of classically chaotic mapstil. In fact, the time-dependent 
Schrodinger equation of these systems is precisely given by Eq.(Q) below. We expect, however, that in the large- N 
limit the results reported here will be general, applicable to any fully chaotic system. This is so because the asyiaptotic 
behaviour of Circular and Gaussian ensembles coincide. As we will discuss here and in more detail elsewhere£3, some 
differences may however remain, in particular in connexion with the discrete nature of time in the Circular ensembles. 

Following DysonO, we start by representing a system not by its Hamiltonian but by an N x N unitary matrix U 
which determines the time evolution of the system according to the equation 

Mn)) = irmo)) (i) 

where n = 0, 1, 2, . . . represents the time in units of an arbitrary interval of time which we fix to unity. The discrete 
nature of the time evolution encoded in this equation is not essential for the moment, and the time will be rescaled 
to a continuous variable later on. The eigenvalues of U are N complex numbers lying on the unit circle 

U\<p k ) = exp(iO k )\<p k ) k = l,---,N. (2) 

Our main observables of interest are the amplitude of return 

N 



C(n) ee (V(O)^(n)) =^> fc | 2 e"^ (3) 
and the transition amplitude to a different state |x) 



fc=i 



N 

T(u)ee ( X |^(n)) =Y. a ^ mBk > ( 4 ) 

k=l 

with ctfc = (<Pk\ip{0)) an d bk = {fk\x)- Due to the discrete and finite nature of the spectrum these are quasiperiodic 
functions of time. However in the large N limit they will be highly oscillatory functions, and our purpose now is to 
study their statistical properties. A statistical analysis of these amplitudes may be introduced by considering |^(0)) 
as a random vector. However, this is not the way we proceed. More relevant to experiments is to fix the initial 
state ve ctor a nd to consider the matrix U as a random matrix belonging to one of the circular ensembles of random 
matricestHlia: (a) symmetric and unitary for systems having time reversal invariance and rotational symmetry or 
having time reversal invariance and integral spin (/3 = 1); (b) unitary if there is no time reversal symmetry (j3 = 2); 
(c) self-dual unitary quaternion for half- integral spin systems with time reversal invariance (/3 = 4). C(n) and T(n) 
thus become random functions of the random variables a,k, bk and 9k whose distribution we wish to study. The 
state |x) will be considered to be statistically independent with respect to |^(0)), and thus the bk's independent from 
the a,k's. For /3 = 4 each eigenvalue in Eqs.(||) and (||) is in fact doubly degenerated (Kramers degeneracy) and 

\a,k\ 2 represents in this case the sum of the square of the amplitudes in the degenerate subspace la^'l 2 + |a^| 2 (and 

analogously a k bk represents the sum a^b^ + a^b^). 

A basic property of circular (as well as Gaussian) ensembles of random matrices is their Hilbert space rotational 
invariance. There are two basic consequences of this invariance. The first one is the statistical independence of the 
eigenvalues with respect to the components of \ip(0)}: P(a, 9) — Pi (a) P2(9), where a — (oi, • ■ • , <2jv), = (<?i, • • • , 0/v) 
and P represents the joint probability density for the corresponding ensemble. The second consequence is the isotropic 
form of Pi 




iM3) = -=— s V Kl -i • ( 5 ) 



The normalization constant Sp n is the surface of a (f3N — 1) dimensional sphere. The coefficients au are real for 
(3 = 1 and complex for j3 = 2 and j3 = 4. We will not need the exact expression for P2(9) but simply note that P^ is 
invariant under a global phase shift of the 9k- Thus the distribution of C{n) depends only on its modulus. 
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II. THE DISTRIBUTION OF THE RETURN AND TRANSITION AMPLITUDES 



We now show that the probability distribution of the real part of C(n) = C\ +i C2 is, in the large N limit, Gaussian. 
Similar steps can be repeated for the imaginary part. From Eq.(^|) the distribution is defined as 



V{d) =< 5 



N 



\a k \ 2 cos(n9 k ) 



fc=i 



> 



where the brackets denote the ensemble average over the amplitudes and phases. Using the integral representation of 
the delta function, the average over the amplitudes reduces to 



N 

Y[ / d 2 a k P 1 (a k )cxp [-i£,\a k \ 2 cos(n0 fc )] 



where P\{a k ) is the distribution of a single component. Because asymptotically the latter distribution is Gaussian, 
then 



I exp 

Z7T 



2 i £ cos(n9 k ) 
N 



> 



where the remaining ensemble average is over the phases. In the large N limit we can expand the logarithm and keep 
only the 1/7V term. The latter equation then reduces to 



nci) « < / f- exp i d d - i Yl cos ^) 

L V fc=i 



> 



i.e. the distribution of the function Ci(n) = J2 k cos(n9 k )/N. This is a function that depends only on the phases (the 
powers of the traces of U). The distribution of the powers of the traces of unitary operators was discussed in Ref. 
|l9| , where it was shown to be Gaussian for low powers. More generally, it has been shownEll that asymptotically the 
distribution of any "linear" function of the form ^2 k f(0k) is Gaussian. Since C\{n) is precisely of that form, this 
then proves that the real and imaginary part of the return amplitude are Gaussian distributed. 

For the transition amplitude the proof is slightly different. For simplicity we limit here to the case (3 = 1 (real 
coefficients); the other symmetries can be treated in the same way. The distribution of the real part of Tin) = T\ +i T2 
is, by definition (cf Eq.(|)) 



P(Ti) =< S 



N 



fe=l 



> 



The use of the integral representation for the delta function and the computation of the integrals over the a's and 6's 
now leads to 



P(Ti 



< 



2tt 



e l « Tl exp 



iVmfl + i^M 



> 



where again the remaining average is over the phases. Keeping the leading order term in 1 /N in the expansion of the 
logarithm it follows that 



< 



2tt 



exp 



N 



2N 2 



k=i 



> 



Finally, computing the integral over £ we arrive to 



P(Ti 



< 



T 2 



2na% i (9) 



: exp 



2o 2 Ti (9) i 
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where 



4 1 (^) = ^E cos2 ( n ^) ■ 

fe=l 

Because in the large N limit the variance cr 2 ^ (0) can be approximated by 

1 f 27T 1 
j o 

it then follows that the distribution of T± is Gaussian with a variance . Likewise, for the imaginary part of T we 
also get a Gaussian distribution with a variance 

1 [ 27T 1 

ov ss / sin (ndi.)d 9 = (1 — S r , n) . 

T2 2-kN J V ' 2N K '° J 

We thus recover, as it should, Eq.(||) below for the variance of T obtained from a direct computation of the average 
return probability. 



III. THE RETURN AND TRANSITION PROBABILITIES 



A. Their average value 

Having established that the asymptotic distribution of the real and imaginary part of the return and transition 
amplitudes are Gaussian distributed, we now consider their first and second moment of those distributions. Of 
particular interest is the second moment of the return amplitude. But before that, let's consider the first moments. 

Since the ensemble average of |afe| 2 is given by < | a.^ | 2 >= J \ak\ 2 Pi(a)da — 1/N, then the average of the return 
amplitude Eq.(||) is 



< C(n) >= — f d6e ind 



Ri(6) 



where R\{9) — N/2n is the level density. Therefore on average the vector \ip(n)) decorrelates immediately from the 
initial state. This is obviously due to the randomization of the phase of \i/)(n)). For the transition amplitude we find 
< T(n) >= 0, Vn. 

The second moment of the distribution of the return amplitude is the probability of return 

N N 

|C(n)| 2 =^K.| 4 + ^| a ,Ha fc |V"^-^. (6) 
fc=i j^k 

As mentioned before, this quantity was considered in earlier work by different authors. For the circular ensembles of 
RMT we obtain 

<lC{n)l >_ W+W) ' ( ) 

which is an exact expression. The function bp\n\ = ^ J Q W d6 cos(n8)T 2 (9) is the Fourier transform of the two-level 



cluster function, the so-called spectral form factorEJ. This function, which for n = is equal to one while for large 
times tends to zero, is responsible for the "correlation hole" mentioned in the introduction. In fact, aside from the 
delta function at the origin, the numerator in Eq.(0) goes from 2/j3 for short times up to its large-time stationary 
value ((3 + 2)j (3. Due to the periodicity of the circular ensembles and the discreteness of time, this result differs from 
the one obtained for the Gaussian ensembles (no convolution present, cf Refs. |Io| , [l6| ). 

Previous works have considered the large N limit of the average return probability by employing the well-known 
large- iV behaviour of the form factor. In order to have an expression valid for any dimension, we have computed the 
form factor for any N, with the result 
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2 ( ^ )/2 / P 2 



A^ ^— ' V p 2 — n 2 

p=l/2 V ^ 



6f ) (r) = (l-r)e(l-r) , (8) 

where we have rescaled the time in units of the Heisenberg time r = n/tjj = n/N, where tjj — 2ir/ A and A = 2ir/N 
is the average level spacing (for (3 = 4 this is the nondegenerate average level spacing). The rescaled time tends to 
a continuous variable in the large N limit. (0(x) is the step function). In a different context, related results where 
obtained in Ref. It is remarkable that for the unitary ensemble the exact and asymptotic form factors coincide. For 
the symplectic case, it is well known that the strong spectral correlations produces a strong " revival" at the Heisenberg 
time (r = 1). Asymptotically this revival is mathematically associated to a divergence of the form factor, since in the 

large- limit it takes the form 64 (t) = lirri/v^oc b^\r) = [1 — r/2 + rln(|l — r|)/4]0(2 — r). However, this asymptotic 
expression does not allow to estimate the exact magnitude of the revival, i.e. if < |C(r)| 2 > is hnite at r = 1 for 
(3 — 4 as N tends to infinity, because both the numerator and the denominator diverge in Eq.fffl). From Eqs.(|^) we 
obtain that at large but finite N the average return probability for (3 = 4 behaves like < |C(r = 1)| 2 >~ In N/N , i.e. 
it is large in a scale 1/N but remains however negligible (no "macroscopic" effect). 

The second moment of the transition amplitude (i.e., the average transition probability < |T(r)| 2 >) differs from 
the average return probability in that it is time independent (no transient regime). In fact, assuming that the state 
|x) is statistically independent from \ip(0)} we obtain for (3=1 and j3 = 2 

<|T(r)| 2 >-<|( X |V(T))| 2 >=l/iV; 0=1,2. (9) 

Analogously we find of course that the average transition probability is 1/(2N) for (3 = 4. 

Before closing this subsection, and generalizing the discussion of Ref. [19] concerning the powers of traces of unitary 
matrices, we would like to include a random walk interpretation of the solutions of the time-dependent Schrodinger 
equation, Eqs.(^-Q), and in particular of the average return probability, Eq.(^). 

If the phases (#1, . . . , On) are considered as random independent variables, their multiples n9k are also random 
independent. Under this assumption concerning the phases, the complex numbers C(n) and T(n) in Eqs.(|^-Q) precisely 
define the position of a particle in the complex plane undergoing a random walk, the length of each of the N steps 
of the walk also being a random variable with a distribution defined by Eq. (^) . In particular, the variance or typical 
square distance traveled by the particle after A" steps will be, for uncorrelated phases 

, ]r ( N5 nfl + {(3 + 2)/(3 

<|C " (n)l >= (N + 2/P) ' (10) 

(the fact that the variance doesn't grow like A" as in the usual random walk process is related to the normalization 
in Eq.(|^) which implies < |afc| 2 >= 1/N, i.e. an average step length which decreases with the number of steps. 
Putting < I ctfc 1 2 >= 1 multiplies the last equation by A^ 2 , thus recovering the usual behaviour). In eq.([lO|) we obtain 
a time-independent variance (since it doesn't depend on n aside the delta dependence at the origin) whose value 
coincides with the asymptotic value of Eq.(0). 

The solutions of the time-dependent Schrodinger equation obtained in the context of RMT may also be thought of 
as a random walk in the complex plane, with the important difference that now there exist correlations between the 
angles of the different steps of the walk. The comparison of Eqs.(^) and ([lO]) indicate that for times shorter than the 
Heisenberg time the variance of this correlated motion is smaller than the variance of the uncorrelated process. This 
is easy to understand intuitively. Indeed, it is well known, and this information is of course contained in the form 
factor, that there exist long range correlations between the (9'st3. For example, the two-point correlation function for 
(3 = 2 exhibits for sufficiently large separations between phases an average repulsion which goes like the inverse of the 
square of the distance. In the random walk process, there exist therefore an average repulsion between the direction 
9k of the kth step and all the other ones. For n = 1, the orientation in the complex plane of one step increases the 
probability for the subsequent steps to go in the opposite direction. This mechanism of course diminishes the variance 
or typical square distance traveled by the particle. 

As n increases, the multiplication of the phases by n and their periodicity contributes to progressively destroy this 
effect, since the repulsion may be transform into attraction by the multiplication and periodization process (depending 
on the relative position of the two phases in the unit circle). At the Heisenberg time, two phases separated by a mean 
level spacing coincide when multiplied by n = N. Then, for times bigger or comparable to the Heisenberg time the 
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product nOk start to behave as a random phase, uncorrelated with respect to the other ones, and the variance of an 
uncorrelated process should be recovered. This is exactly what is observed in Eq.(Q) for (3 = 1 and (3 — 2, since the 
function bp(n) monotonically reaches the value zero at t « 1. For (3 = 4, however, the level repulsion is so strong 
that it produces a very rigid structure for short distances: the probability to observe equidistant phases 

9 k = 2-Kk/N (11) 

for k = 1, 2, 3, . . . is considerable enhanced. This strong local order, which is practically absent for (3 = 1 and (3 = 2, 
decreases with k and becomes irrelevant for k = 5 — 10 (see for instance Fig. 6 in Ref. ||). In the random walk process, 
all steps which approximately satisfied Eq.(|ll]) are in phase since n0 k = N9 k = 2nk. This produces the observed 
behaviour of the variance for (3 — 4 at the Heisenberg time, which grows like In N/N. 

Finally, we comment on the absence of transient regime in the average transition probability (Eq.(^|)). The difference 
between Eqs.(^|) and (Q) is that in the former |dfc| 2 is a real random number, while a k bk is a complex one (or real with 
a random sign for (3 = 1). This means that the prefactor in Eq.(Q) contains a random phase on it, which of course 
suffices to destroy all the correlations between the (9's. The transition amplitude in quantum mechanics may then be 
thought of as a standard random walk process with uncorrelated phases. The variance obtained (Eq.(|)) coincides, 
aside from the normalization, with the usual result (multiplying by -/V 2 to recover the usual normalization we get a 
variance which grows linearly with the number of steps). 



B. Their distribution 



The distribution of the return and transition probabilities is obtained straightforwardly from that of the amplitudes. 
For example, introducing the rescaled variable -q = \C(t)\ 2 /(2 < \C(t)\ 2 >), it follows from the Gaussian distribution 
of C(t) that the distribution of the variable i] is given by 

V{rj) = cxp(-7y) . 

An analogous equation is valid for the transition probability. This distribution was in fact conjectured by Prang 
in his study of the statistical properties of the form factor. As was discussed there, the Poissonian distribution of 
rj implies that |C(r)| 2 and |T(t)| 2 are non self-averaging quantities, since the mean and the variance of r\ are both 
equal to one and therefore the fluctuations compare to the mean do not tend to zero in the large- N limit. In order 
to observe this quantities, some additional smoothing (over time or over the parameters of the system) is necessary. 
In fact, Eq.(^) represents an interferent sum of oscillating terms. The shortest period in the sum corresponds to an 
eigenphase difference of order 7r, giving T m in = O(l) ~ t erg (see the concluding remarks for a discussion concerning 
terg)- On the contrary, the longest period is associated to differences in the eigenphases of the order of the mean level 
spacing, T max ~ 2ir/ A = tn ■ In order to keep the time-dependent structure of the curve, a reasonable smoothing is 
over an interval 5t such that t erg <C St -C tjj, thus reducing the fluctuations by a factor of order 1 / \5t. Averages 
over the initiaLstate or parameters of the system may also be considered, like it was done in the experimental studies 

of < |c(r)i 2 >ma. 



IV. CONCLUDING REMARKS 



In Eq.([y), one iteration of the unitary matrix corresponds to the time it takes to the system to explore (without 
filling) the available phase space. This is because the state U\ip(0)) can roughly be at any point of Hilbert space (U 
is an arbitrary rotation in that space). This time, called the ergodic time t erg , can be relatively large in a particular 
physical system as compared to other time scales involved in the problem, meaning that the system is not able to 
immediately explore the available phase space. For example, in disordered metallic systems a diffusive motion takes 
place before the particle reaches the boundary of the sample. If the typical sample size is L and the diffusion constant 
D, then t erg cx L 2 / D (Thouless time) which may be much larger than the elastic mean free_time. The diffusive motion 
is known to produce non- universal deviations from the RMT in the form factor for t < t erg c3. These non- universalities 

can also manifest themselves as a system-dependent transient in Eq. (^|) . Non- universal features will appear in b^ (t) 
for short times in general for any dynamical systemcj, where t erg can roughly be identified with the period of the 
shortest periodic orbit. As extensively studiedl 2 !, in a scale of the order of t erg the function b^\r) may show a 
series of peaks at multiples of the shortest periodic orbits whose amplitude decreases as exp(— Xt), where A is the 
classical Lyapunov exponent. For longer times, it follows from the Gaussian distribution of the return amplitude with 
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a variance of order 1/N that statistically the probability of revivals (i.e. a return probability of order one) in chaotic 
systems is exponentially small. 

Although a natural choice (relevant for experiments) of 1^(0)) is a wave packet, it is important to note that the 
relaxation process described here is not restricted to such a class of initial states. Our results are general, valid for 
any initial state (except of course the eigenstates). In particular, \ip(0)) could be a state completely delocalized in the 
basis we are using to study the motion. 

In conclusion, we have studied the relaxation process in the time evolution of complex quantum systems. Aside 
from the short time non-universal features already mentioned, the relaxation process described here is expected to 
be universal, applicable to a large class of physical systems including those whose classical limit is chaotic, disor- 
dered metallic mesoscopic samples, high lying excitations of many body systems, and equivalent problems in optics or 
acoustics. The 1/N scaling factor of the fluctuations amplitude is a purely normalization factor, proportional to the 
typical component of the initial state |?/>(0)) on the chaotic eigenstates of the system. For example, for maps defined 
on a compact phase space, N is the dimension of the Hilbert space, while for chaotic billiards N is the volume of the 
cavity. Rescaling the time by the Heisenberg time we obtain a universal description, expected to be independent of 
any specific feature of the system like for example the distribution or density of impurities or the Lyapunov exponents. 
In Ref. ^ numerical results for the transition probability for several chaotic maps were reported. These were found 
to be independent of the system considered and they all nicely agree with the present results. 

We would like to thank Y. Alhassid, O. Bohigas, O. Legrand, O. Martin and N. Pavloff for fruitful discussions and 
comments. We also thank the Institute for Nuclear Theory at the University of Washington for its hospitality and 
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